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ABSTRACT 

We investigate the kinetics of constant-kernel aggregation which is aug- 
mented by either: (a) evaporation of monomers from finite- mass clus- 
ters, or (b) continuous cluster growth - i. e., condensation. The rate 
equations for these two processes are analyzed using both exact and 
asymptotic methods. In aggregation-evaporation, if the evaporation is 
mass conserving, %. e., the monomers which evaporate remain in the sys- 
tem and continue to be reactive, the competition between evaporation 
and aggregation leads to several asymptotic outcomes. For weak evapo- 
ration, the kinetics is similar to that of aggregation with no evaporation, 
while equilibrium is quickly reached in the opposite case. At a critical 
evaporation rate, the cluster mass distribution decays as /c -5 / 2 , where 
k is the mass, while the typical cluster mass grows with time as t 2 / 3 . In 
aggregation-condensation, we consider the process with a growth rate for 
clusters of mass k, which is: (i) independent of k, (ii) proportional 
to k, and (iii) proportional to k^, with < \x < 1. In the first case, 
the mass distribution attains a conventional scaling form, but with the 
typical cluster mass growing as tint. When oc k, the typical mass 
grows exponentially in time, while the mass distribution again scales. 
In the intermediate case of oc /c M , scaling generally applies, with the 
typical mass growing as $1/(1-/0. We also give an exact solution for the 
linear growth model, Lk oc k, in one dimension. 



P. A. C. S. Numbers: 02.50.-r, 05.40.+j, 82.20.-w, 82.20.Mj 
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I. INTRODUCTION 



Aggregation, fragmentation, and condensation underlie a variety of non-equilibrium 
phenomena in nature [1-6] . In systems where only one of these mechanisms is operative, the 
time-dependent cluster mass distribution generally evolves to a scaling form in which the 
basic variable is the ratio of the cluster mass to the typical mass. These scaling forms have 
been obtained by exact solutions, numerical simulations, and by direct consistency checks 
of the scaling description. By these efforts a general understanding has been developed 
for the connection between microscopic reaction details and macroscopic features of the 
cluster distribution. 

These approaches have also been successfully applied to processes where the mech- 
anisms of aggregation, fragmentation, and condensation are simultaneously active. One 
important example is aggregation in combination with fragmentation, a process which 
arises naturally in reversible polymerization [7]. Since the basic elements of aggregation 
and fragmentation are manifestly opposed, their combined effect generally leads to an 
equilibrium state in a closed system. Characterizing the approach to and the detailed na- 
ture of this equilibrium has been of basic interest [7-11]. Detailed balance considerations 
can generally be applied to determine the nature of the equilibrium state. If Kij denotes 
the aggregation rate of Ci + cj — > Ci+j, where Ck is the concentration of clusters of mass 
k, and Fij denotes the fragmentation rate of Ci+j C{ + Cj, then detailed balance gives 
KijCiCj = FijCi + j where Ck is the steady-state concentration of /c-mers. 

On the other hand, the situation where the fragmentation matrix has mostly zero 
elements naturally arises in polymer chain growth kinetics [12,13]. For example, /c-mers 
may be unstable to break-off of monomers ("evaporation"), while all other fragmentation 
events are forbidden. By the nature of these restrictions, the process may be viewed 
as aggregation in a "dry" environment. For this system, the aforementioned detailed 
balance argument can no longer determine the equilibrium state. One of our goals is 
to investigate the kinetics of this composite aggregation-evaporation process in the rate 
equation approximation. We generally consider the situation where the aggregation rate 
is independent of the masses of the two incident clusters. We further restrict ourselves 
to the interesting case of mass conservation, where the monomers which have evaporated 
remain in the system and continue to participate in the reaction. If the evaporation is 
sufficiently strong, its dominance over the effects of aggregation results in an equilibrium 
whose properties are determined analytically. In the opposite case where aggregation 
dominates, the typical cluster mass increases indefinitely. However, at a critical value of 
the evaporation rate, the mass distribution decays as /c _5//2 , while the typical cluster mass 
grows with time as t 2 / 3 . These intriguing features emerge from both exact solutions to the 
rate equations and asymptotic arguments. 

In a complementary direction, we also explore the kinetics of a combined aggregation- 
condensation process in which each cluster grows at a specified mass-dependent rate, in 
addition to the aggregation. This underlying growth can be viewed as arising from a 
uniform flux of monomers which permeate the system. In this sense, the composite reaction 
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can be viewed as aggregation in a "wet" environment. Such a process naturally arises in 
various contexts, such as the growth of breath figures and thin film deposition and growth 
[6,14,15]. For breath figures, in particular, theoretical models have generally considered 
the growth rate of each droplet to be a specified function of time. This ultimately gives 
rise to much faster droplet growth in the composite aggregation-condensation process. In 
contrast, we consider a droplet growth rate which depends only on the cluster size, and 
not explicitly on the time. We also employ the rate equations, a description which may 
not necessarily apply to the stationary droplets which are characteristic of breath figure 
systems. Within our approximation, we find a number of surprising features in the time 
evolution of the cluster mass distribution for a mass-independent aggregation rate. When 
the condensation rate is independent of the cluster mass, the distribution of cluster masses 
obeys conventional scaling, but with the typical mass growing as t Int. On the other hand, 
for a condensation rate which is proportional to the mass, the typical mass grows as e*. For 
a general mass-dependent condensation rate which is proportional to k^ (with < \i < 1), 
the typical mass grows as 

£V(i-m). For all three cases, that is, < fi < 1, scaling generally 
applies and this approach further predicts that the mass distribution decays exponentially 
with mass. 

A detailed treatment of these two cases of aggregation-evaporation and aggregation- 
condensation is presented in the following two sections. 



II. KINETICS OF AGGREGATION-EVAPORATION 



The rate equations for mass conserving aggregation-fragmentation are, 

(la) 

+ [L k+1 c k+1 (t) - L k c k (t)} + S k ,i y^LjCj(t). 



Here the overdot denotes the time derivative and c k (t) is the concentration of clusters of 
mass k at time t. In this equation, the first two terms account for the gain and loss of k- 
mers due to aggregation, respectively. The prime on the sum in the gain term indicates the 
restriction of mass conservation, i + j = k. In the evaporation process, a k-mer produces 
a (k — l)-mer and a monomer at a rate L k = A/c M . The gain and loss of /c-mers because of 
evaporation are described by the third and fourth terms of Eq. (1(a)). Finally, the last term 
accounts for monomer production as a result of evaporation. (If monomers were removed 
from the system by evaporation, the last term would be absent. This non-conservative 
process has trivial kinetics in which the mass disappears exponentially in time.) 

Let us now determine the conditions for which the system either reaches equilibrium 
or evolves ad infinitum. We focus on the case where the aggregation rates ify are all equal 
(constant-kernel aggregation) and where the evaporation rate is independent of the cluster 
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mass (fj, = 0). Physically, this latter rate is appropriate for linear polymers with evapora- 
tion possible only at the chain ends. Although one can easily generalize the discussion to 
mass-dependent aggregation and evaporation rates, their relative influences are typically of 
different orders of magnitude. Consequently, it is relatively straightforward to anticipate 
whether an equilibrium or a scaling distribution obtains. However, for mass-independent 
aggregation and evaporation, the competition between these two influences is delicate and 
gives rise to surprisingly rich kinetic behavior. 

With the assumptions of constant reaction kernel and size-independent evaporation, 
the rate equations simplify to 

1 / 

d k(t) = 2^ K c ^) c iW " ^K Cl (t)+\ (c k+1 (t) - c k (t)) + A 4,i Ci(t). (lb) 

i,j i=l i=l 

The evaporation rate A and reaction rate K can be absorbed by redefining the concentra- 
tions and time by c k {t) — > 2Xck(t) and t — > t/XK, leading to 

oo oo 

ck(t) = J2' a(t)cj(t) -2c k (t)^2a(t) + (c k+1 (t) - c k (t)) + d^J^^it)- ( lc ) 

i,j i=l i=l 

For simplicity, consider a monomer-only initial condition, c k (t = 0) = A48 k ,i- In this case, 
the total initial mass is the only control parameter, with a large mass corresponding to a 
small evaporation rate and vice versa. 

To gain insight into the kinetics, it is helpful to first write the equations for moments 
of the mass distribution, M n (t) = ^ k> ik n c k (t). By straightforward manipulations on 
Eq. (1(c)), these moments obey the equations, 

M (t) = -Mo(t) 2 + Mo(t)- Cl (t), 
Mi(t) =0, 

M 2 (t) =2(M 1 2 -M 1 + M (t)), 

M 3 (t) = 3Mi + 3(2Mx - 1)M 2 (£), ( 2 ) 
M 4 (t) = 2M (t) - 4Mi + 6(M 2 (£) + M 2 {tf) + 4(2M 1 - l)M 3 (t), 
M 5 (t) = 5Mi - 10M 2 (t) + 10M 3 (t) + 20M 2 (t)M 3 (t) + 5(2M 1 - l)M 4 (t), 



For the monomer-only initial condition, M n (t = 0) = A4 for all n. From the equations for 
M n (t) for n > 2, it is clear that the typical mass and higher moments grow indefinitely, if 
the initial mass is sufficiently large. In the complementary case, however, an equilibrium 
state is possible. 

More complete information about the kinetics can be obtained by analyzing the rate 
equations themselves. For this purpose, we introduce the generating function 

oo 

g(z,t) = J2ck(t)z k . (3) 
k=i 
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By multiplying the rate equation for each Cfc(t) by z k and summing over all k, the generating 
function obeys 



g(z,t)=g(z,t) 2 -2g(z,t)M (t) + [ - Cl {t) ) - g(z,t) + M (t)z. (4a) 



Here we use the equivalence g(z = l,t) = Mo(t), with Mo(t) the cluster number density. 
As is often the case in these types of systems, it is more convenient to consider a modified 
generating function, h(z, t) = g(z, t) — g(l, t) = g(z, t) — M (t), in which the value at z = 1 
is subtracted off. This generating function satisfies 

h(z, t) = h 2 (z, t) + ^ h(z, t) + (1 ~^ 2 M (t). (46) 

z z 

While we have been unable to solve this differential equation in general, the time 
independent solution is 

Kz) = ^ (l - y/1 - AM z) , (5) 

with h(z) = h(z, t = oo), Mq = Mo(t = oo), and the sign of the radical is fixed by requiring 
that h(z) — > — Mq as z — > 0. Once the value of M is specified, the equilibrium solution, 
as well as the conditions for equilibrium to exist can be determined. From Eq. (5), we 
conclude that equilibrium exists whenever 4M < 1, while if 4M > 1, the power series 
representation of the generating function diverges and equilibrium is not reached. The 
unknown quantity M can be related to the initial mass by the requirement that z^| z=1 , 
which is the total mass of the system, equals the initial mass Ai. This leads to the condition 
M = (1 - VI - 4M ) /2, or equivalently, M = M(l - M). Since M must be in the 
range (0, 1/4) for equilibrium to occur and since Mq must be an increasing function of Ai, 
the corresponding constraint on the initial mass is M. G (0, 1/2). 

The equilibrium properties of the cluster size distribution can be obtained by expand- 
ing the generating function in Eq. (5) in powers of z for 4M < 1. Thus the equilibrium 
cluster mass distribution, cjt = Ck(t = oo), is given by 



1 

Cfc 



4v^F 



(6) 



where T(z) is the gamma function. When 4M < 1 (equivalently, M < 1/2), the asymp- 
totic behavior of Ck is dominated by (4Mo) fc and the mass distribution decays exponentially 
in k. On the other hand, when 4Mq = 1 (A4 = 1/2), the mass distribution has the power- 
law form 

3 T(k — 2) , _ 5 /2 /P7 n 

Cfc = — f= — oc k 1 . 7) 

8^ T(A; + 2) v ; 

The behavior of the moments of the mass distribution reflects the above two possi- 
bilities. Since all positive moments approach equilibrium values for TV! < 1/2, a recursive 
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solution of the moment equations, Eq. (2), gives 



M 
Mi 



M(l-M), 
M, 



M 2 



M 



M 3 



1-2M' 
3M(1-M) 
2(1 -2M) 3 



M(l + M) 
2(1 -2M) ' 



(8) 



On the other hand, for the limiting case of M = 1/2 (M = 1/4), the power-law form of 
the cluster mass distribution leads to finite values of the moments M n for n < 3/2 and 
diverging values for n > 3/2. 

When M. > 1/2, an equilibrium solution of the rate equations does not exist and 
the transient behavior is of basic interest. To understand this behavior, it is helpful to 
first consider the equation of motion for M 2 (t). For M 2 (£) to remain positive, M (t) must 
approach a constant whose value is greater than M(l — M) as t — > oo (Eq. (2)). This 
observation suggests that the cluster mass distribution breaks up into two components for 
M. > 1/2. In the scaling component, the typical cluster mass M 2 (£) diverges as t — > oo. 
Conversely, the remaining non-scaling (equilibrium) component arises from the continued 
production of monomers by evaporation which attempts to maintain the number of clusters 
Mo and the mass M. at their equilibrium values of 1/4 and 1/2, respectively. Thus the 
scaling part of the mass distribution should contain the remaining mass of M. — 1/2, while 
the number of clusters in this sub-population should decay to zero. 

Since the non-scaling component of the mass distribution is given by Eq. (7), the 
corresponding contribution to M n (t) is given by M n (t) ~ J2k<t^ n ~ 5 ^ 2 ~ t n_3 / 2 for n > 
3/2, while M n (t) — > const. < oo for n < 3/2. These estimates are based on the observation 
that the mass distribution approaches the steady-state form Ck ~ /c _5//2 behind a leading 
edge which is growing linearly in time, while there are essentially no clusters of larger mass. 
To determine the relative contribution of the scaling component of the mass distribution to 
the moments, notice, from Eq. (2), that the typical cluster size M 2 (t) oc t. This suggests 
that the ratio x = k/t is the appropriate scaled mass with which one can express the 
scaling component of the mass distribution as 



The prefactor in Eq. (9) is needed to have a finite mass in this scaling part of the dis- 
tribution. If we are interested only in the asymptotic behavior, then from Eq. (9), the 
contribution of this component of the mass distribution to the n th moment is 



c k (t) =t- 2 $(x). 



(9) 




where 




(10) 
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As a result of these estimates, it is evident that the moments with n > 1 are dominated 
by the scaling component of the mass distribution. We therefore restrict ourselves to n > 1 
and insert the above scaling expression into the moment equations. Again, retaining only 
the asymptotically relevant terms, Eqs. (2) are transformed into a system of equations for 
the amplitudes A n 

n-2 , s 

(n-l)A n = (2M-l)nA n - 1 + J2[ n i ) A i A n-i, n>3. (11) 

i=2 ^ ' 

To determine A2, we use the third of Eqs. (2) together with the expected asymptotic value 
Mo = 1/4. This gives A2 = 2(A4 — i) 2 . To solve for the remaining A n we introduce the 
generating function 

3=2 J ' 

It is then straightforward to transform the system (11) into the differential equation 

dA ( , . , 1 x 2 



Solving Eq. (13) yields 



and therefore 



dz y^M--}. (13) 



A (14) 



M n (t) = „< [ M - i) f~\ (15) 



From the moments given by (15), it is evident that the scaling component of the cluster 
mass distribution is 

<*(') = jM^^[-JJ^i\- < 16 ) 

We therefore arrive at the same expression for the mass distribution as that which arises 
in pure aggregation with a constant reaction rate [1-4]. The only difference is that the 
mass which comprises the scaling component is equal to M. — \. 

To summarize, in the weak aggregation regime, M. < ^, the mass distribution ap- 
proaches the equilibrium form of Eq. (6) at an exponential rate in time. In the complemen- 
tary strong aggregation regime, M. > -^i the typical mass grows linearly in time. However, 
there is an anomalous enhancement in the small-mass tail of the mass distribution which 
is of the form given in Eq. (7). This residue arises from the continued re-introduction of 
monomers into the system by evaporation. 

At the critical point M. = -|, a novel temporal behavior can be anticipated in which 
the approach to equilibrium occurs at a power-law, rather than an exponential rate. Let 
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us therefore hypothesize that the number of clusters, M (t), approaches its equilibrium 
value of 1/4 as t~ v . Employing this assumption in the moment equations gives the series 
of relations M 2 (t) ~ t x ~ v ', M 3 (t) = ±4^ _ tj M 4 (t) ~ t 3 " 2iy , M 5 (£) ~ t 3 "^, etc.. Since 
the exponents of successive reduced moments should be equidistant within a scaling de- 
scription, the condition M4/M3 ~ M 3 /M 2 , e. implies f = |. This suggests the general 
formula, 

M^)-^" 1 (17) 

for n > |, while for n < | the moments approach finite values. 

One can refine this analysis to obtain the amplitudes as well. Writing the expected 
asymptotic behavior as M n (t) = B n t^~ x for n > | and substituting into the moment 
equations, the B n satisfy recursion relations for integer n. Defining Bq via Mo(t) — | = 
B t~ 2 / 3 , we then find B 2 = 6S , B 3 = §, and 



/2w 



2 ' 

n-2 



7t — Z /■ x 



(18) 



for n> A. While we could not solve Eq. (18) generally, we can obtain partial information 
by a generating function approach. Indeed, introducing 



00 _ 



*(*) = £ if**" 1 ( 19 ) 



i=2 J ' 

transforms Eqs. (18) into a Riccati equation 

d 4- = B 2 + B z~ 2 / 3 + -. (20) 
dz 4 

While Eq. (20) is insoluble when B 2 ^ 0, several characteristics can be established without 
an explicit solution. In particular, the structure of Eq. (20) suggests that a solution 
exists for z < z c and it has form of a simple pole, B(z) ~ (z c — z)~ x , in the vicinity of 
z = z c . Inversion of this singularity yields the large-n behavior for every third amplitude, 
B 3n — 2(3n)!z~ 2n , while -63^+1 and -E?3 n +2 are asymptotically negligible compared to B 3n . 
The location of the singularity, z c , depends on the unknown amplitude So and should be 
established from a more comprehensive treatment which, in turn, is equivalent to finding 
a complete solution of Eq. (20). 

III. KINETICS OF AGGREGATION-CONDENSATION 

(a) Mass Independent Growth Rate (L^ = const.) 

We now consider the complementary reaction where bimolecular aggregation (with 
mass-independent aggregation rates) , is supplemented by single-cluster growth in which a 
cluster of mass k grows at a rate L k oc We wish to understand how this additional 
growth influences the kinetics of the underlying aggregation within the rate equations. Let 
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us first investigate this composite aggregation-condensation process for a size-independent 
growth rate, L k = A. The rate equations become 



oo 

Ck(t) = J2' ^(t) Cj (t) - 2c k (t) ^q(£) + A (c fc _i(t) - c k (t)) . (21) 

i,j i=l 

To gain qualitative insight into the asymptotic behavior, we begin by solving for the first 
few moments of the mass distribution. We then present a complete solution for the mass 
distribution, from which the asymptotic behavior may be extracted. 
From the above rate equations, the moments evolve according to 

M (t) = -M (t) 2 , 
Mi(t) = AMo(t), 

M 2 (t) = 2M 1 (t) 2 + X(2M 1 (t) + M (t)), (22) 



subject to the initial condition M n (t = 0) = 1 for all n. We also set the initial density 
equal to unity so the condensation rate A is the only control parameter. Solving for the 
moments successively yields 

M (t) 



i + r 

Mi(t) = l + Aln(l + £), 

M 2 (t) = l + Aln(l+t) + (4A-2A 2 )(l + t)ln(l + t) (23) 
+ 2A 2 (1 + t) ln 2 (l + t) + 2(1 - A + A 2 )t, 



Although the exact expressions for M n (t) become cumbersome as the index n grows, the 
asymptotic behavior is simply 



M n (t)^ n \X n t n - 1 (\nt) n . 



(24) 



To solve the full rate equations we again introduce the generating function g(z, t) = 
J2T=i c k(t)z k , which reduces an infinite set of rate equations (21) to the differential equa- 
tion 



g(z,t) = g{z,t) 2 - 



g(z,t) + X(z-l)g(z,t). 



(25) 



1 + t 

Notice that gi^z.t) -1 satisfies a linear inhomogeneous differential equation whose solution 
is 

ze \(z-l)t r f t dr 



1-z 



-, -1 



,A(*-l)i 



o (1 + r)' 



(26) 
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We now determine the asymptotic behavior associated with this generating function. This 
information will also suggest the asymptotic form for the mass distribution for general 
rates of aggregation and condensation, a system for which an analytical solution cannot 
be found. First, notice that the densities Ck(t) which make a non-zero contribution to 
the generating function are those whose mass is in the range < (1 — z)k < 1. If we 
further assume that (At) -1 , then the integral in Eq. (26) can be approximated 

by replacing the exponential by unity; this is asymptotically correct over the domain of 
integration. The generating function now becomes 



z e 



\tz 



g (z,t)^(l + t)- 2 e- xt ^— . (27) 

1 — z 

Expanding Eq. (27) gives, 

Ck+i(t) = ^-t -2 e -A *, for < k < At. (28) 

Thus for the range k <C At, the mass distribution is Poissonian; however, the distribution 
cannot be written in the conventional scaling form t~ a ^(k/t^). 

On the other hand, for sufficiently large k, the mass distribution does exhibit scaling. 
To determine an appropriate mass scale we expand the exponent e A ( 2-1 - )T , compute the 
integral in the right-hand side of Eq. (26), and then asymptotically balance the various 
terms. The nontrivial scaling limit arises when (1 — z) ~ (tint) -1 and suggests that the 
appropriate scaling variable is 

(= (l-z)Xtlnt, (29) 

instead of the original variables z and t. 

In the scaling limit, t — > oo and 1 — z — > with £ kept fixed, the generating function 
simplifies to 

g(z,t)^t- 1 (l + C)- 1 , (30) 
and the mass distribution approaches the scaling form 

c k (t) ~ <Kt)Q(x), with x=j^- v (31) 

with the prefactor <j)(t) and scaling function $(x) to be determined. Making use of Eqs. (29) 
and (31), we express the generating function in terms of <f>(t) and as, 

~ Xtlnt (p(t) / <$>{x)e~ xC dx. 
Jo 

Finally, by comparing Eqs. (30) and (32) and performing the inverse Laplace transform, 
the prefactor <f>(t) and the scaling function $(x) are 

^ ) = At4it' and ®^ =e ~ X - ( 33 ) 
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The scaling solution of Eqs. (31) and (33) agrees with asymptotic expression for the mo- 
ments (24). 

For completeness, we also investigate the large-mass tail of the mass distribution, 
k ^> Xt In t. The analysis is similar to that given above so we merely cite the result 

^-{i-jL)'- (34) 

Thus the mass distribution function does not scale in both the small- and large-mass tails. 
Formally, the condensation process governs the small-mass tail of the distribution, as well 
as the overall mass. Conversely, the form of the distribution in the scaling region and in 
the large-mass tail is determined solely by the aggregation process. 

(b) Growth Rate Proportional to the Mass (L^ oc k) 

We now consider the case of a condensation rate which is linear in the mass, i. e., the 
rate of Cf- — > Ck+i equals Xk. The corresponding rate equations are 

oo 

Ck(t) = Yl' c ^) c jW " 2c fc (t) ^c,(t) + X[(k - l)c fc _i (t) - kc k (t)] . (35) 

i,j i=l 

Employing the generating function g(z,t), Eq. (35) becomes 

-g(z, t) = g(z, t) 2 - 2g(z, t)M (t) + Xz(z - 1) —g(z, t). (36) 

Notice that the number of clusters M (t) = g(z = l,t) still satisfies Mo(t) = — M (t) 2 ; 
hence, Mo(t) = (1 + Introducing again the modified generating function, h(z,t) = 

g(z,t) — M (t), transforms Eq. (36) into a linear equation for /i(z,t) _1 

-h(z, t)- 1 + Xz(l - z)—h{z, t)~ 1 + l = 0. (37) 

By introducing w = jliij^, Eq. (37) becomes a first-order wave equation in the 
variables (w,t). This equation further simplifies by transforming from (w,t) to u = t + w 
and v = t — w to yield 

Ak(«,«)-' = 4 (38) 

with solution /i(w,w) _1 = — ^ + f(v). Here f(v) is fixed by the initial conditions. For 
example, for monodisperse monomer-only initial conditions, we obtain 

h( u , v)- 1 = -1 - ^ - e~ Xv = -l-t- e~ Xt -^. (39) 

By now expanding g(z, t) = h(z, t) + in powers of z the exact concentrations are, 

-\t / e ~\t N k ~ 1 



Ck[t) = (TTip V 1 " T+7, 1 (40) 
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In the scaling region, k — > oo and t — > oo with a; = k/te xt finite, Eq. (40) has the asymptotic 
form 

Cfc(t) ~t- 2 e- At e- x . (41) 

Notice that the scaling solution for the mass distribution in both the cases of ji = 
and fx = 1 may be written as 

, , M (t) 2 _ . , M (t) , 

Although the mass density, Mi(t), has a very different time dependence for these two 
cases, the respective scaling functions are the same and in fact identical to that in pure 
aggregation with a constant aggregation rate [1-4] . 

(c) General Mass Dependent Growth Rate (L;t oc 

We now turn to aggregation-condensation for a general homogeneous mass-dependent 
cluster growth rate, = A/c M , with < /z < 1. Although it does not appear possible 
to solve the governing rate equations, one expects the scaling form of Eq. (42) to hold in 
the scaling region. This assumption reduces the problem to finding the first two moments, 
Mo(t) and M\(t). The former task is trivial since the condensation process does not alter 
the evolution of the number of clusters, and Mo(t) = (1 + t)~ . On the other hand, 
Mi(t) is determined by M\{t) = AM M (t) which is coupled to an undetermined moment. 
Fortunately, in the long-time limit we can use the scaling form (42) to estimate M M (t) as 

(43) 

Thus asymptotically M x {t) ~ AT(1 + ^)t /i_1 Mi(t)^, which may be solved to yield 

Mi(t) with A=[\(l- n)Y(n)}^ . (44) 

It is instructive to compare this asymptotic result for the typical cluster size, 

with a naive estimate provided by considering growing, but noninteracting, i. e., non- 
aggregating, clusters. This latter estimate follows from S(t) = A>S(t) M , which implies 

S(t) = A t^, with A = [A(l -/z)] 7 ^ 7 • (46) 

Therefore the system with growing but non-interacting droplets provides the correct expo- 
nent of the time dependence for the droplet growth rate in the interacting system. However, 
the corresponding prefactor A is slightly smaller than that of the interacting system. 
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(d) Scaling Approach for Low Spatial Dimension 

For diffusion-controlled aggregation, the above mean-field approaches are typically not 
applicable for spatial dimension d < 2 (see, e. g., [16,17]). However, for the aggregation- 
condensation process with a homogeneous growth rate, oc /c M (0 < \i < 1), it is possible 
to infer partial results for d < 2 by applying scaling and exploiting known results. In 
particular, for diffusion-controlled aggregation with a mass-independent cluster diffusivity, 
the density of clusters (which is not altered by the condensation process) is [16,17] 

( f~ d / 2 d < 2- 
M.W~{' ¥> ' HI (47) 

We now again assume that asymptotically the mass distribution approaches the scaling 
form 

Cfc (i)~W $(l) with x =M^ k , (48) 
v ; M x {t) v ; Mi(t) ' v ; 

with a general (i-dependent scaling function ^^(x). The mass density, Mi(t), is determined 
from Mi(t) = AM M (t), where the moment M M (£) is estimated to be (following the steps of 
the previous subsection), 

poo 

M/j,(t) ~ Mo(t) 1_M Mi(t) M / dxx^ d (x). (49) 

Jo 

Ignoring numerical factors we solve for the mass density to obtain 

{t^h~i, d<2, < fj, < 1; 
t^lnt, d = 2, 0<//<l; (50) 
ln 2 t, d = 2, fi = 0. 

Finally, combining Eqs. (47), (48), and (50) yields 



t !-m * <& d (x), d < 2, < /i < 1; 
c fc(*) ~ W"t^ hit $ 2 (a0, d=2, 0<//<l; 
1 t" 2 $ 2 (^), d= 2, // = 0. 



(51) 



in all three cases, the scaling variable is x = /ct -1 /^ 1- ^. 

(e) Exact Solution in One Dimension for Lk oc k 

We now solve the aggregation-condensation process in one dimension when the clus- 
ter growth rate is linear in the mass (Lk oc k) and the diffusivity of each cluster is mass 
independent, a case particularly amenable to exact analysis. The model is defined on the 
real line which is populated by point clusters which: (i) diffuse with a mass independent 
diffusivity D, (ii) grow with a rate proportional to their masses, and (iii) aggregate irre- 
versibly whenever two clusters meet. A convenient way to treat this problem analytically 
is to introduce [18] the quantity Pk(x,t), which is the probability that the total mass of 
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all clusters contained in an interval of length x at time t is equal to k. The evolution of 
the Pk(x,t) are uncoupled and each such function obeys the diffusion equation [18]. If we 
account for condensation, the Pk(x,t) now evolve according to 

d 

t) = 2D—P k (x, t) + X [(k - l)P fc _i - kP k ] . (52) 

The second term clearly accounts for the change in mass in the interval due to the cluster 
growth process. The cluster concentrations, c k (t), are then found from 

c k (t) = ^-P k (x = 0,t). (53) 

Note also the sum rule, Ylk^o ^k{ x i = 1> which allows us to focus on the P k s with k > 1. 
To solve this system, we introduce the generating function, 

oo 

V(z,x,t) = ^z k P k (x,t), (54) 



k=l 



to transform Eq. (52) into 



dV d 2 V dV 

- = W — + Mz -,)-. (55) 



Note that the obvious boundary condition, P k (x = 0) = Ofco, translates to the condition 
on the generating function, V(x = 0) = 0. If initially the system is composed of monomers 
which have a Poissonian distribution of positions with density unity, the corresponding 
initial condition is 

k 

P k (t = 0) = ^-e- x , V(t = 0) = e --( 1 --) - e ~ x . (56a) 
k\ 

Although the variable x is defined only on the half-line x > 0, it proves useful to consider 
Eq. (55) for — oo < x < oo, and impose the initial condition 

V(t = 0) =e x -e x{1 ~ z \ for x < 0. (566) 

With this antisymmetric initial data, the boundary condition V(x = 0) = is manifestly 
satisfied. 

Employing the auxiliary variables u = t + j In and v = t — ^ In introduced 
earlier, we further simplify Eq. (55) to the diffusion equation 

^= D ^ (57) 
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subject to the initial conditions 



V\ = e - x /( 1+e ~ X ^-e- x , 

I u= — v 

V\ = e x - e x ^ 1+e ~ Xv \ 

I u= — v 



for x > 0, 
for x < 0. 



Solving Eq. (57) with the above initial conditions yields 



V(u, v, x) = 



— f 

Hu + v) Jo 



^4nD(u + v) 



exp 



x 



exp - 



1 + e Xv 
AD{u + v) 



- exp - 



(^ + 2 

AD{u + v) 



Returning to the original variables we find 

1 



V(z,x,t) = 



V8nDi Jo 



x 



exp 



exp 



1 + e \tl^z 



8Dt 



exp 



8Dt 



Combining Eq. (60) with Eq. (53) gives 



exp 



1 

1 + e At(i^) 



- 1 



(58) 



(59) 



(60) 



(61) 



Eq. (61) is the formal solution to the problem - for example, the concentrations can be 
found by expanding the right-hand side in powers of z. However, explicit expressions for 
the concentrations Ck(t) are cumbersome, even for small k. In contrast, the moments of 
the mass distribution have simpler form. For example, 



Mo(t) 



y/32ir(Dt) 



Wf <^ exp (-^) (1 - e ~ £) 



1 



V2nDt' 



(62) 



in agreement with the known result [19]. Furthermore, the mass density is given by M\(t) = 
e xt . This is expected since M\{t) satisfies M\{t) = \Mi{t). These two moments provide a 
useful consistency check of our exact results. 



IV. SUMMARY AND DISCUSSION 



Our primary result was to elucidate the broad range of phenomenology which arises 
from the combined effects of: (a) aggregation with evaporation, and (b) aggregation with 
condensation. In the former case, the interesting situation is that of mass independent 
rates of aggregation and evaporation in a mass conserving system. An equilibrium state 
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is reached for sufficiently strong evaporation, while the kinetics is essentially identical to 
that of pure aggregation when the evaporation is relatively weak. At a critical evaporation 
rate, there is power-law kinetics in which the typical cluster size grows as t 2 / 3 , while the 
mass distribution decays with mass k as /c -5 / 2 . An essential ingredient in these results is 
that the aggregation and evaporation rates are of the same order, so that their competition 
leads to interesting manifestations. Related transition behavior was obtained for combined 
aggregation-evaporation by Virgil et al. [10], but with both a mass-dependent aggregation 
rate (proportional to the product of the cluster masses) and evaporation rate (proportional 
to the cluster mass) . For this latter system, the effects of evaporation and aggregation are, 
in some sense, of the same order of magnitude, leading to a transition between equilibrium 
and gelation for a critical value of the ratio of the two rates and also distinct power-law be- 
havior at the transition. Another not entirely unrelated example where microscopic effects 
influence the approach to equilibrium occurs in the one-dimensional reversible reaction 
A + A <-> A [20] . In this case, the relaxation has different functional forms depending on 
the ratio of the initial density to the final equilibrium density. 

For aggregation with condensation, the cluster growth enhances the aggregation, as 
expected. When the growth rate is independent of cluster mass, this enhancement is 
relatively weak, however, and the typical cluster mass grows with time as tint, compared 
to a linear growth for aggregation with no condensation. Conversely when the growth 
rate is proportional to the mass, the typical mass grows exponentially in time. In the 
intermediate case of a growth rate for clusters of mass k given by oc with < \i < 1, a 
scaling approach indicates that the typical cluster mass grows as t l ^ l ~^\ while the scaling 
function in the mass distribution is a pure exponential for < \i < 1. 

Our aggregation-condensation model is related to breath figures [14,15] as well as to 
other models of droplet growth and coalescence [6]. In these latter systems, the growth 
rate of individual droplets is generally a specified function of time, rather than of the 
mass. If one assumes that the mass of individual droplets grow t a , then the typical cluster 
mass in the combined aggregation-condensation system grows at a faster rate of t@, with 
(3 = Da/(D — d). Here D is the spatial dimension of the droplets and d (which must be 
smaller than D) is the dimensionality of the substrate. It would interesting to understand 
what relation, if any, exists between these results and our model of aggregation with explicit 
mass-dependent cluster growth. 

The rate equation approach is expected to provide the correct asymptotic behavior for 
aggregation-condensation when the spatial dimension d is greater than 2 [16,17]. To extend 
our understanding to d < 2, a simple-minded approach was developed which suggests that 
conventional scaling continues to apply. Furthermore, we derived an exact solution for 
aggregation-condensation in one dimension with a growth rate linear in the cluster mass, 
Lk oc /c, and with a mass-independent diffusivity. The scaling predictions are in satisfying 
agreement with this exact solution. 
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